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Nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) techniques have 
supplied us with unique information on local incoherent phenomenon as well as low frequency mag- 
netic response of electronic systems. Some of recent advances in the studies of high-Tc cuprate 
superconductors are based on site-selective measurements of the local electronic state and the lo- 
cal electron spin dynamics. High-Tc superconductors YBa2Cu307_5 (the optimal Tc = 93 K) and 
YBa2Cu408 (Tc = 82 K) with double-layer Cu02 planes have been studied by NMR and NQR 
techniques. From the intensive studies at the early stage, two-dimensional antiferromagnetic spin 
fluctuations and pseudo spin-gap in the magnetic excitation spectrum turned out to be key ingredi- 
ents to characterize the novel electronic states, associated with the occurrence of superconductivity. 
Wipeout effect is a new key to characterize superconductor-insulator boundary in the lightly doped 
regime and impurity substitution effects. What happens near in-plane impurities and inside vortex 
cores in such novel electronic states? Nonmagnetic impurity Zn decreases Tc of the d-wave super- 
conductivity because of the breakdown of Anderson theorem. From the careful NMR and NQR 
studies, it was concluded that the pseudo spin-gap in the normal state is robust for the substitution 
of nonmagnetic impurity Zn and magnetic impurity Ni. In the vicinity of a superconductor-insulator 
boundary induced by impurity, however, the pseudo spin-gap is suppressed. Collapse of low energy 
pseudo spin-gap due to the heavily substituted impurity Zn was also evidenced by neutron scattering 
measurements. For the impurity-substituted samples and in the vortex state of pure samples, the 
characteristic NMR and NQR spectra were observed. The local enhancement of magnetic correla- 
tion near the impurity Zn and inside the vortex core was observed from the frequency distribution 
of nuclear spin-lattice relation time T\. Scanning tunneling spectroscopy detects local density of 
states of electrons near Zn in the superconducting state. NMR and NQR experiments detect the 
local magnetism near Zn and inside the vortex cores. 

PACS numbers: 



I. INTRODUCTION 

Nuclear magnetic resonance (NMR) and nuclear 
quadrupole resonance (NQR) are powerful techniques 
to characterize local incoherent phenomenon in materi- 
als [E whereas neutron and X-ray scattering tech- 
niques are sensitive to coherent phenomenon Q. Mi- 
croscopic studies using the NMR and NQR techniques 
have provided us with rich information on crystalline im- 
perfection, dilute exchange spin networks, Friedel oscil- 
lation, Ruderman-Kittel-Kasuya-Yosida (RKKY) oscilla- 
tion, and Kondo screening effect in materials [1, 0, d, 0| • 
Some of recent advances in the NMR and NQR stud- 
ies of high-Tc cuprate superconductors are based on site- 
selective measurements of the local electronic state and 
the local electron spin dynamics. 

YBa2Cu30j; (6.4< y <7) (Y123y) system is a widely 
known high-Tc superconductor. Y123y has two inequiv- 
alent crystallographic Cu sites; bilayer plane-site Cu(2) 
and single chain-site Cu(l). Figure [T] shows the crystal 
structures of the Y123 family. The maximum supercon- 
ducting transition temperature Tc is 93 K at the optimal 
oxygen concentration of y = 6.92~6.96. For oxygen defi- 
ciency of 6.4< y <6.6, the Tc ~ 60 K phase is known to 
show a typical underdoped electronic state with pseudo 
gap and pseudo spin-gap. Figure [5] is a typical phase 
diagram of the Y123 |8|. 

Non-stoichiometry in compounds gives rise to inhomo- 
geneous broadening in NMR and NQR spectra, which of- 
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FIG. 1: Crystal structures of pure undoped YBa2Cu306.o, 
optimally doped superconductor YBa2Cu307, and naturally 
underdoped superconductor YBa2Cu40g. 



ten prevents us from extracting intrinsic electronic prop- 
erties. YBa2Cu408 (Y124) with bilayer Cu02 planes and 
double CuO chains has Tc — 82 K. It is a stoichiomet- 
ric and naturally underdoped compound. The oxygen 
concentration is hard to be removed. Because of chemi- 
cal stability and of lack of oxygen disorder effects, since 
the first observations of Cu NQR NMR and NQR 

studies have been intensively carried out for Y124 as well 
as the Y123 family to reveal intrinsic electronic proper- 
ties 



2 



500 



400 



300 



200 



100 



' 1 ' 1 


1 1 1 1 1 

YBa^CujOy 


antifeiTomagnetic 1 


- 


1 . 1 


superconducting 
1,1, 



6.0 6.2 6.4 6.6 6.8 7.0 
Oxygen concentration y 



FIG. 2: Phase diagram of Y123y reproduced from Rossat- 
Mignod et al. @|. Tn and Tc are plotted against oxygen con- 
centration y. 



The aim of this article is to present some of the re- 
cent advances of NMR and NQR studies of high-Tc su- 
perconductors Y123 and Y124. The article is orga- 
nized to show (II) the present understanding of magnetic 
phase diagram with respect to doped hole concentration, 
(III) impurity-substitution effects, (IV) impurity-induced 
NMR and NQR relaxation, and (V) vortex core state in 
the superconducting mixed state. Recent scanning tun- 
neling microscopy (STM) and spectroscopy (STS) stud- 
ies of high-Tc superconductors have supplied us with vital 
information of local density of states of electron ifl^ [Tst- 
Local microscopic magnetism near impurity and vortex 
core, which we are concerned about, is the matter that 
NMR and NQR can clarify actually. 



II. MAGNETIC PHASE DIAGRAM 

A. TWO DIMENSIONAL RENORMALIZED 
CLASSICAL REGIME 

The parent undoped compound YBa2Cu306 is a Mott 
insulator and bilayer-coupled Heisenberg antiferromag- 
net with a Neel temperature Tn = 415 K. The plane- 
site Cu(2)^+ carries an unpaired moment S = 1/2. The 
chain-site Cu(l)^"'" has a closed shell. The paramagnetic 
state is in two dimensional renormalized classical regime, 
where the antiferromagnetic correlation length £ diverges 
exponentially with cooling toward T = K [3, EE Ell 



C/a « 0.5exp(27rpjT), 



(1) 



where ps is a spin stiffness constant proportional to an 
in-plane exchange coupling constant J. The value of J 
is estimated to be about 150 meV il!l|. Because of a 



spin S' = 1/2 system, there is no single-ion anisotropy 
due to electrostatic crystalline potential. The spin wave 
dispersion in the Neel state is nearly gapless EmI; which 
leads to high Cu nuclear spin-lattice relaxation rate at 
low temperatures below Tn [1^. Thus, the Tn higher 
than room temperature can be considered to result from 
the strong in-plane J with XY anisotropy and three 
dimensional interactions of J±_ as /cbTn « J±{^/a)'^ E3- 
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FIG. 3: (a) High temperature Cu nuclear spin-lattice re- 
laxation rate l/Ti of parent insulator YBa2Cu306 reproduced 
from Pozzi et al. [S^l and of naturally undersoped YBa2Cu40g 
reproduced from Machi et al. [23| . (b) The nuclear spin-lattice 
relaxation rate l/Ti calculated by SCR (self-consistent renor- 
malization) theory from weakly antiferromagnetic to nearly 
antiferromagnet in two dimension reproduced from Sli]. To is 
a spin fluctuation energy. T/To is a reduced temperature. 



The renormalized classical behavior in Cu nucle ar sp in- 
lattice relaxation was first observed by Imai et al. |22l | for 
La2Cu04. The measurement of high temperature spin 
dynamics is a highlight in the successful NMR spin-echo 
experiments. Subsequently, such a renormalized classical 
behavior in NMR relaxation was observed by Pozzi et al. 
for Y1236 [23i], Cao.ssSro.isCuOa (Tn = 539 K) [21, by 
Fujiyama et al. for PrBa2Cu408 (Tn = 230 K) [p" and 
by Thurber et al. for three- leg ladder Sr2Cu305 [2g |. 

Figure [3] (a) shows the high temperature ^^Cu nuclear 
spin-lattice relaxation rate 1/Ti of YBa2Cu306 repro- 
duced from Pozzi et al. [23| and of naturall y und erdoped 
YBa2Cu408 reproduced from Machi et al. |l32l |. 
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B. SUPERCONDUCTOR-INSULATOR 
BOUNDARY 
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FIG. 4: Magnetic phase diagrams of Yi_a;Caa;Ba2Cu30j; re- 
produced from Singer and Imai [35l |. 

As the hole concentration increases, Tn is reduced 
and finally the superconductivity emerges. Accord- 
ing to a field theory applied to a non-linear sigma 
model [itI . [H, [2^ and to a modified spin wave theory (30l| 
in the insulating regime, a spin frustration effect destroys 
the coherency of an antiferromagnetic network and finally 
a paramagnetic state or disordered gapped state should 
emerge through a two dimensional quantum critical point 
(QCP). According to the SCR (self-consistent renormal- 
ized) theory for a two dimensional spin fluctuation model 
in the itinerant regime, the static staggered magnetisa- 
tion in a weakly antiferromagnetic state is reduced and 
finally a nearly antiferromagnetic state emerges through 
an itinerant two dimensional QCP [31]. 

Figure [3] (b) shows the nuclear spin-lattice relaxation 
rate calculated by SCR (self-consistent renormalization) 
theory from weakly antiferromagnetic to nearly antiferro- 
magnet in two dimension reproduced from (3l| . Although 
the development of the spin dynamics with hole doping 
from Y1236 to Y1248 in Fig. [3] (a) is similar to that via 
the QCP in Fig. [3] (b), the actual cuprates change from 
insulator to superconductors through a Mott transition. 
What happens at the boundary from insulating side to 
superconducting side? 

Wipeout effect on NMR is a key to character- 
ize the physics at the boundary of the Y123 fam- 
ily as well as La2-a;Srj;Cu04 [33, fss!, (33|. Sudden 
disappearance of Cu NMR and NQR signals is ob- 
served near the superconductor-to-insulator boundary of 
Yi_a;Caa;Ba2Cu30jy. Figure|4]shows a magnetic phase di- 
agram of Yi_2:Ca2:Ba2Cu30y [3^. In general, the wipe- 
out is recognized as an ill-defined physical quantity. Nev- 
ertheless, it well characterizes what happens near the 



phase boundary. An abrupt wipeout effect is different 
from conventional slowing down phenomenon, where the 
nuclear spin relaxation should become fast. The wipeout 
effect without any slowing down effect on the observable 
NMR/NQR signals leads us to conclude that the electron 
spin fluctuation spectrum should consist of two compo- 
nents in frequency space. One is a main spectrum, which 
may be called a fast mode. The other is the induced low 
frequency spectrum, similar to a central mode in struc- 
tural phase transition, which may be called a slow mode. 
In Fig. m the existence of the slow mode is evidenced by 
ligand Y NMR relaxation measurements [s^ . Because 
of a weak ^^Y nuclear-electron coupling constant [s^, the 
fast nuclear relaxation component being unobservable by 
Cu nuclear spin-echo is observed by Y NMR. 
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FIG. 5: Two components of Cu electron spin fluctuation 
spectrum. The slow mode is inferred from wipeout effect on 
NMR and NQR. Such a mode near superconductor-insulator 
boundary might be associated with Zhang-Rice bound sin- 
glets. 

The spin fluctuation spectrum with two energy scales 
reminds us of a phenomenological itinerant-localized dual 
model [37I [38j . The actual slow mode is inhomogeneous 
in space and in frequency. In most cases, Mott transition 
is first order transition. Thus, such a slow mode may be 
associated with the character of first order transition. 

Holes are doped into the plane-site oxygens [39| . In the 
superconducting regime, the doped holes do not have in- 
dependent spin degree of freedom, which is indicated by a 
single scaling relation of "'^^0(2, 3) with ^•^Cu(2) Knight 
shifts [401 . The leading spin component of conduction 
band is single [H, HO, S^l . The doped holes do not carry 
isolated free moments but are considered to form Zhang- 
Rice bound singlets of 3(PL [4l|. The Zhang- Rice sin- 
glets are localized in the semi-conducting regime at low 
temperatures. The phenomenological slow mode may be 
associated with the localized Zhang-Rice singlets, which 
cause the wipeout effect on NMR/NQR. 

Can we observe NMR/NQR signals associated with 
the Zhange-Rice singlets in any cases? One may find a 
few examples of high frequency Cu NQR signals in some 
doped cuprates. For the hole-doped ladder and linear 



4 



chain oxides of Sri4Cu2404i '43'| and Cao.85Cu02 [44| . 
Cu NMR/NQR signals near the Zhang-Rice singlets are 
observed in the chain sites. In spite of the singlets, the 
nuclear spin-lattice relaxation time is relatively short. 
For Li-substituted La2Cu04 ^45] and NaCu02 ^4^], non- 
mganetic Cu^+ NQR signals are observed to be in low 
spin electronic states. The observed magnetism of Cu^+'' 
of 3(PL is different from Cu'^+ so far. One should proceed 
further site-selected NMR studies to observe the local 
electronic state of Zhang-Rice singlets in doped oxides. 



C. NORMAL-STATE PSEUDO SPIN-GAP 

Normal-state pseudo gap and pseudo spin-gap are 
characteristic anomalies of the underdoped superconduc- 
tors. The existence of normal state gap is inferred from 
the above-Tc decrease of nuclear spin-lattice relaxation 
rate divided by temperature 1/TiT [H, \M, \M, H \M, 
Isoj . 1/TiT is the square of hyperfine coupling constant 
times the wave vector averaged low frequency dy- 
namical spin susceptibility x" (q, t^) (wave vector q and 
frequency uj) (51| . Inelastic neutron scattering confirms 
the low frequency suppression of x" (q, ^) for underdoped 
Y123 [11]. The effect of the pseudo spin-gap is also rec- 
ognized in transport properties [s^ . 
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FIG. 6: Pseudo spin-gap of Y(Bai_a;Laa;)2Cu30^ reproduced 
from Matsumura et al. [s^l and Auler et al. [5J|. The hole 
concentration P^h per Cu02 sheet is defined by Cu^"*"^"''. 

Figure [H] shows the typical hole doping de- 
pendence of pseudo spin-gap temperature Tg for 
Y(Bai_a;Laa;)2Cu30y, which is defined by the maximum 
temperature of Cu jH, The effect of the 

opening of pseudo spin-gap is not a phase transition but 
a crossover phenomena. The Ts decreases as the hole 
concentration per Cu increases and then coincides with 
Tc at slightly overdoped concentration. Although the ex- 
istence of overdoped regime of Y123y is known for single 
crystals, it was unclear for polycrystalline samples. Ya- 
mamoto et al. found the method of a solid state reaction 
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FIG. 7: Hole doping dependence of pseudo spin-gap tem- 
perature Ts as a siglet-RVB formation temperature Trvb in 
two dimensional t-J model factorized by spinon and holon 
operators reproduced from Tanamoto et al. ifjfi] (a), and as 
an onset temperature (a mean-field T^^^) of d-wave supercon- 
ducting fiuctuations reproduced from Yanase and Yamda [s^l 
(b). 



and heat treatment to synthesize the overdoped polycrys- 
talhne Y1237 with Tc = 88 K [l^. 

One should note that the above- Tg behavior of 
is of Curie- Weiss type but not Korringa 
law [sl, [5^1. Although the non-Korringa behavior in 
the overdoped regime is not widely recognized, the 
Curie- Weiss type l/^^TiT is also confirmed for over- 
doped Lai.8oSo.2oCu04 [5^, Lai.77So.23Cu04 [58], and 
HgBa2Cu04+5 [59|. The high temperature electronic 
state in the overdoped regime is not conventional 
Landau-Fermi liquid. There is no phase boundary in 
the high temperature electronic state at the optimally 
doped-to-overdoped regimes. 

Since the pseudo spin-gap effect was not obvious in 
single layer La2_a;Sr2Cu04, it had been suspected that 
such a gap is not essential for high Tc. Does the 
pseudo spin-gap result from a spin singlet paring in the 
adjacent Cu02 layers as a property peculiar to dou- 
ble layer compounds [60j or from an inherent property 
in a single Cu02 plane? It had been a long term 
problem. Distinct conclusion was obtained from sin- 
gle layer HgBa2Cu04+5 [6l|, The pseudo spin-gap 
is an inherent nature of the high-Tc superconductors, 
irrespectively of the number of Cu02 planes per unit 
cell. The doping dependence of Tg of the single layer 
HgBa2Cu04+5 is different from that of the double layer 
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HgBa2CaCu206+5 [6^. The Tg as a function of the hole 
concentration depends on the number of Cu02 planes per 
unit cell and on the shape of underlying electronic Fermi 
surface. 

Two theoretical explanations have been proposed. In 
the two dimensional t-J model with spinon-holon de- 
composition technique, the pseudo spin-gap tempera- 
ture Tg is regarded as the onset temperature Trvb 
of singlet RVB (resonating valence bond) formation of 
spinous [s^, E^, 'GdI . The real Tc is given by Bose-Einstein 
condensation temperature TbeCj leading to the under- 
doped regime. In the existing calculations, Jrvb is given 
by a second order phase transition temperature. Figure[7| 
(a) shows a numerical calculation of Tryb as a function of 
doped hole concentration [6^. The doping dependence 
of Trvb depends on the contour of a basal Fermi sur- 
face. For each high-Tc family with different Fermi sur- 
face, Trvb shows a different doping dependence. 

In the two dimensional superconducting fluctuation 
theory, the pseudo spin-gap temperature Tg is regarded as 
the onset of enhancement of d-wave superconducting fluc- 
tuations and a mean- field Tc. The actual Tc is reduced 
by the strong superconducting fluctuations so that the 
underdoped regime appears. Thus, the mean- field Tc is a 
crossover temperature. Fi gur e [71(b) shows the mean-field 
Tc and the suppressed Tc [67[ 



III. IMPURITY-SUBSTITUTION EFFECT 

A. PAIR BREAKING AND SUPERUNITARY 

Magnetic and nonmagnetic impurity substitution ef- 
fects have been studied for the high-Tc superconduc- 
tors, where the impurity is substituted for the plane- 
site Cu(2). The impurity Zn (spinless S — Q) induces 
Curie magnetism, a local staggered spin susceptibility in 
the normal state, an in-gap antiferromagetic dynamical 
spin susceptibility, and a virtual bound electronic state 
around Zn in the superconducting state. 

In strong coupling superconductors, the actual su- 
per conducting transition temperature Tc is determined 
by competition between pairing effect and depairing ef- 
fect [68]. In the phonon-mediated superconductors, the 
depairing effect due to thermal phonon competes the par- 
ing effect due to virtual phonon. For s-wave supercon- 
ductors, non-magneitc impurities do not suppress super- 
conductivity because of Anderson's theorem, unless they 
induce Anderson localization effect. For d-wave super- 
conductors, non-magneitc impurities suppress supercon- 
ductivity. This is a natural consequence of breakdown 
of Anderson's theorem (toI . Fzij . The superconducting 
state of the high-Tc superconductors is characterized by 
dj.2_y2--wave pairing [69|. 

The criterion of scattering strength of impurities is 



given by the value of 7„ 



7„ = TriVpu, 
7,„ = ttNfJS, 



(2) 



where Np is the density of states of electrons at the Fermi 
level, a nonmagnetic potential depth u, and an exchange 
coupling constant J of a magnetic impurity S with con- 
duction electrons. A weak scatterer can be treated by 
lowest order Born approximation as |7n.m| < 1. A strong 
scatter can be treated by unitarity limit as |7n.m| > 1. 

In an s-wave superconductor, weak magnetic scatterers 
suppress superconductivity: 



ATc = 0.5Tr^Ci^pNFJ^SiS+l)/kB, 



(3) 



where Cimp is in-plane impurity concentration. This is the 
Abrikosov-Gorkov's formula. In a d-wave superconduc- 
tor, weak magnetic scatterers suppress superconductivity 
by 



ATc = 0.257r2chnpiVF J^S'(S' + 1)/A:e 



(4) 



In the lowest order Born approximation, the c?-wave su- 
perconductor is twice robust than the s-wave supercon- 
ductor for the destructive magnetic impurities [73 . [t^ . 

For the magnetic impurity effect in s-wave supercon- 
ductors, one can neglect nonmagnetic potential scatter- 
ing with each magnetic impurity, because of Anderson's 
theorem. However, the potential scattering also affects 
d-wave superconductivity as well as the magnetic impu- 
rity scatte ring . Thus one has to take into account both 
effects [7l[74[,[zl- 

Taking into account the impurity scatterings of clas- 
sical spins (Sz = iS) within a i-matrix approximation, 
Ohashi derived the theoretical expression of Tc (cimp <Cl) 



a dr>. 



-wave superconductor [76], 



ATc 



1 



4fcBiVF 



1 



+ (7n-7m)^ l+(7n+7m)^ 



(5) 

For 7m = 0, eq.(l5|) gives ATc in the case of nonmagnetic 
impurity Zn [77, 78, 79|, [13|, whereas 7„ = gives ATc 
in the case of magnetic impurity without nonmagnetic 
potential scattering. In the limit of — > and 7„ — > 0, 
eq. ([5]) leads to the expression of ATc in the lowest order 
Born approximation |72l | . In addition, tv™ — > cxd or 7„ — > 
oo corresponds to the unitarity limit [8l[: 



ATc = 



*--imp 



(6) 



If cancellation of 7m 
for Ni, one obtains 



in the unitarity limit occurs 



ATc - 



'-'imp 
SfcBiVF ' 



(7) 



which is two times smaller than that due to pure non- 
magnetic scattering. 
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FIG. 8: Tc versus impurity concentration z for various high- 
Tc superconductors with M=Ni or Zn. The sohd (open) sym- 
bols are the Ni(Zn) doping dependence of Tc- The data are 
adopted from l^.[8a ] for YBa2(Cui-zMz)307 (upward trian- 
gles), [8l,|8i,l86[ for YBa2(Cui-zMz)40s (circles), and 
for Lai.85Sro.i5Cui-zMz04 (downward triangles). The de- 
crease of Tc by Ni doping for Y123 is most slow among these 
compounds. Ni doping dependence of Tc for "as-synthesized" 
or "quenched" Y123, adopted from |82|. The solid curves are 
guide for the eyes. The decrease of Tc by Ni doping is larger 
in "quenched" Y123 than in "as-synthesized" one. 

In Fig. [HI impurity substitution effects on Tc are plot- 
ted against the concentration z for various high-Tc su- 
perconductors with impurity M=Ni or Zn. The solid 
(open) symbols are the Ni(Zn) doping dependence of Tc. 
The data are adopted from [S^, HS] for optimally doped 
YBa2(Cui_^M^)307 (upward triangles), [H, IH, jsll 
for naturally underdoped YBa2(Cui_2M^)408 (circles), 
and [131 for Lai.85Sro.i5Cui_zM204 (downward trian- 
gles) . The degree of suppression of Tc by Zn and Ni impu- 
rities depends on the hole doping level in Y123y [sl]. For 
less concentration of oxygen, ATc/ Acimp becomes larger. 

The decrease of Tc by Ni substitution for the optimally 
doped Y1237 is most slow among these compounds. Us- 
ing a reduced oxygen partial pressure technique, Adachi 
e< al. succeeded in controlling Tc of YBa2Cu3_2;Nia;07_5 
with optimal oxygen concentration [s^ . They synthe- 
sized two series of YBa2Cu3_a;Nia;07_5 samples with the 
different Tc's per Ni concentration. We call the sam- 
ples with Tc denoted by black upward triangles "as- 
synthesized'' ones, because they were synthesized in flow- 
ing oxygen gas without quenching treatment, and the 
other samples with lower Tc denoted by red upward 
triangles "guenc/iecf' ones, because they were the as- 
synthesized samples once again fired and quenched in a 
reduced oxygen atmosphere at 800 °C. In Fig. [51 Ni sub- 
stitution dependence of Tc is plotted for "as-synthesized" 
or "quenched" Y123 [13]. The decrease of Tc by Ni 
is larger in "quenched" Y123 than in "as-synthesized" 
one. Ni prefers the higher coordination of oxygen atoms. 
The plane-site Cu(2) is located in the pyramid with five 



oxygen ions, whereas the chain-site Cu(l) is coordinated 
with two, three, or four nearest neighbor oxygen ions. 
The two series of samples with different Tc suggest that 
the distribution of Ni-substituted sites over Cu(l) and 
Cu(2) sites is changed through synthesis under the re- 
duced oxygen atmosphere. 

From the Cu(l) and Cu(2) NQR measurements, it 
turned out that the heat treatment in reduced oxygen 
atmosphere results in a redistribution of the Ni ions in 
YBa2Cu3-^Ni:r07-5 The plane-site ^^Cn{2) nu- 

clear spin-lattice relaxation for the quenched sample was 
faster than that for the as-synthesized sample, in con- 
trast to the ^■^Cu(l) relaxation that was faster for the 
as-synthesized sample. This indicates that the density of 
plane-site Ni(2) is higher in the quenched samples, con- 
trary to the chain-site Ni(l) density which is lower in 
the quenched samples. The Ni substitution for the chain 
Cu(l) site is a reason why the suppression of Tc by Ni 
impurity is so slow in the optimally doped Y1237. In 
passing, Co impurity is also substituted both for Cu(l) 
and Cu(2) sites in the 1237, which was verified by Cu(l) 
and Cu(2) NQR by Kohori et al ^ prior to the Ni 
study m- 

Monthoux and Pines had an insight into the impu- 
rity substitution effect on the host spin fluctuation spec- 
trum, which yields the pairing interaction of high Tc, as 
well as the life time of superconducting quasi-particles 
of Y123 [131 • They pointed out that the observed sup- 
pression of Tc per Zn is too strong to be ascribed only 
to a simple potential scattering effect on quasi-particle 
lifetime: 

unitary 

« 5(K/%) X (lOOcunp). (8) 

According to the Cu NMR results for Y123 at the 
time [qH ]. Zn impurity induces a large residual Knight 
shift and Korringa-like behavior in the superconducting 
state and non-exponential nuclear spin-lattice relaxation 
in the normal sate, while Ni does not alter the normal 
and the superconducting properties. Taken together with 
these NMR results for Y123, the impurity effect was clas- 
sified into 

J Zn : superunitary, , , 

1 Ni : subunitary. ^ ' 

The substitution of Zn impurity induces the changes in 
the host spin fluctuation spectrum, while Ni does not 
change it but induces the depairing effect on the quasi- 
particle life time. They proposed the strong suppresssion 
of Tc due to a magnetic vacancy correlation function with 
a vacancy operator ft,(r) of an effective radius of ^: 



Xzn(r - r') = xpurc(r - v')h{v)h{v'). (10) 

The magnetic correlation is suppressed over a distance 
of ^. However, this correlation function is not compati- 
ble with all the subsequent NMR and neutron scattering 
results. 
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The fact of suppression of superconductivity by non- 
magnetic impurities has been known from early stages. 
By the Cu NMR techniques, Ishida et al. found that 
Zn impurity in Y123 increases the residual Knight shift 
in the superconducting state [9l[. Although they did 
not observe the Zn-neighbhor Cu NMR signals, they 
found that the observed Cu NMR signals are shifted 
by Zn impurities. From the residual Cu Knight shifts, 
they inferred a quasi-partcile energy band formation 
of Zn-induced virtual bound states. Strong scattering 
impurity potential in the unitarity limit must induce 
a virtual bound state in the d-wave superconducting 
state [H, m, H S S [13, H [9i|. Individual virtual 
bound states around Zn an d Ni imp urities are actually 
observed by STS technique |lO0l . llOll. It rema ined to be 
solved whether the Zn strengthens [l02l Il03| or weak- 
ens [m local magnetic correlation around Zn in the op- 
timally doped and overdoped systems and how the mag- 
netic correlation near and away from Zn is changed or 
unchanged below and above T^. Zn-neighbor NMR and 
NQR measurements were highly desired. 

As for the optimally doped Y1237, a part of Ni im- 
purities is substituted for the chain-site and the amount 
of Ni substitution for the chain site can be controlled 
to some extent [zi, The Cu NQR studies of the 

Ni-redistributed Y1237 as well as Ni-substituted Y1248 
revealed that Ni also induces nonexponential NMR relax- 
ation and suppresses Tc (tgI. I104| . Further, the detailed 
analysis of nonexponential relaxation revealed that no- 
exponential relaxation due to Zn indicates the change of 
the host spin fluctuation spectrum but that due to Ni 
does not necessarily indic ates t he change of t he spin fluc- 
tuation spectrum {ll, M, ^M, [IqS, OM, JM] ■ 



B. SITE-SELECTIVE '^^Y NMR NEAR Zn 
IMPURITY 

The local magnetic correlation is not suppressed near 
Zn impurity. Mahajan et al. succeeded in observing the 
^^Y NMR si gnals near Zn impurities in Zn-substituted 
YBa2Cu30y [llO| |. It is an outstanding site-selective ob- 
servation in the NMR studies of impurity-substituted 
high-Tc superconductors. Figure [5] shows *^Y NMR 
spectrum (inset figure) for YBa2(Cuo.98Zno. 02)306.6 and 
the si te-se l ected Ti near and away from Zn, reproduced 
from [not lllll |. The Zn impurity induces Curie type 
magnetism near the impurity site and enhances the nu- 
clear spin-lattice relaxation rate. Impurity-site NMR 
studies have also been intensively carried out for Li- 



substituted Y123y 



been mtensive ly c arried out 

[MfTTl[Ti3l . [nlfT]l fiT6t. 



INELASTIC NEUTRON SCATTERING FOR 
Zn-SUBSTITUTED Y123 
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FIG. 9: *^Y NMR spectrum (inset) and temperature depen- 
dence of *^Y nuclear spin-lattice relaxation rate divided by 
temperature l/TiT for YBa2(Cui_2Zn^)3 06.64 at T = 100 K 
and at a magnetic field applied alon g ah plane. All the data 
are reproduced from Mahajan et al. [llOl] . 



tions of wave vector q and frequency w. Kakurai et al. 
observed a large change of the dynamic al sp in suscepti- 
bihty of Zn-substituted YBa2Cu3 06.35 |ll7( . Not only 
low energy but also high energy spectrum is changed 
by Zn substituion. Subsequent neutron scattering stud- 
ies revealed for slightly-underdoped and optimally doped 
samples that not a simple shift of the spectrum weight 
but the deve lopm ent of an in-gap state is induced by Zn 
substitution [ll8| |. 
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FIG. 10: Zn-induced change in dynamical spin susceptibility 
X" (q, ijj) around an antiferromagnetic wave vector Q mea- 
sured by inelastic neutron scattering experiments. All the 
data are r eproduced from Kakurai et al.{a) for underdosed 
Y123 [li3l and Sidis et al.{h) for optimally doped Y1237 ful. 



Inelastic neutron scattering experiments allow us to 
see the dynamical spin susceptibility x" (q? '-^) as func- 



The large change at low frequency was confirmed by a 
complete wipeout effect of Cu NQR signals by Kobayashi 
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et al. for underdoped Y123 [119|. A partial wipeout ef- 
fect by Zn impurities has already been obser ved by Yam- 
agata et al. for optimally doped Y1237 |l02{ . One should 
note that the relation between the in-gap state and the 
wipeout effect on Cu NQR for Zn-substitution is parallel 
to that for superconductor-insulator boundary. 



IV. ZERO FIELD CU NQR STUDIES OF 
IMPURITY-INDUCED EFFECTS 



IMPURITY-INDUCED NUCLEAR 
SPIN-LATTICE RELAXATION 
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In contrast to ^^Y NMR, Zn-induced effect on Cu 
NMR and NQR spectra are not apparently s eparable 
near and away from Zn impurity [qJ ll2Cll Il2l( . Then, 
there had been a pessimistic view that no one can ex- 
tract separable information of local magnetism from 
Cu NMR and NQR experiments. However, it turned 
out that the application of an impurity-i nduced nuclear 
spin-lattice relaxation theory jl23l Il24| enables us to 
see such local info r mation near the impurit ies from Cu 
NQR [li, H [Toi [Toi, [loE [l03, Eoi, [103. AIso, the 
finding of Zn-induced Cu NQR signals promotes the Cu 
NQR studies of th e local an tiferromagnetic correlation 
near the impurities [l09l Il22j | . 

Let us briefly explain the impurity - induced nuclear 
spin-lattice relaxation theory jl23l . [124] . This theory has 
been applied for dilute Heisenberg insulators and dilute 
alloys. Assuming only the random distribution of impuri- 
ties, a stretched exponential recovery curve is derived for 
any magnetic insulators and itinerant systems, because 
it is concerned with the kinetics of nuclear moments. 

Figure [TT] is a top view of a Cu02 plane with magnetic 
impurities. Solid red circles are the Cu nuclear sites. The 
relaxation process of a nuclear moment at a i-site consists 
of two parts: 

(fluctuations)host + ^^(fluctuations)imp, 

where the summation is taken over a distribution of im- 
purities in a configuration. If one can observe the re- 
covery curve p{vi, t) of one nuclear spin at a ith site in 
Fig. it must be a single exponential function 



(r,,t) = exp[-t{(;^)host +J2 rrL \ J']' 



exph(^)host]n 



cxp[ 



(11) 



where 1/Ti(r) — A(r)^S'imp(r, u;„) is a nuclear spin- 
lattice relaxation rate occurring by a nuclear-impurity 
coupling. Here, r is a distance between a nuclear and 
an impurity-induced moment, A{r) is a coupling con- 
stant between the nuclear spin and the impurity-induced 
moment, 5'imp('', w„) is the impurity-induced spin-spin 



FIG. 11: A square lattice Cu02 plane with magnetic impuri- 
ties. Solid red circles are Cu nuclear sites. The nuclei within 
a radius (wipeout radius) of of the impurity are unobserv- 
able in conventional resonance frequency and time scale. Blue 
large circles indicate the wipeout regions. 



correlation function [75|, |97|, |98|. a;„/27r is a nuclear 
resonance frequency. For an isolated local moment on 
the impurity site, we have a longitudinal direct dipole 
coupling j4(r)^ cx a nd a two- dimensional RKKY 

interaction A{rf cx l/r^ [l23l. [l2^. For an impurity- 
induced moment with a local A{r) w ^(0), we have 
S'imp(f, w„) cx r^'* and e~^/^ i s a correlation length, 
and a cx is a constant) [lOj|. Instead of including 
all the direct and indirect nuclear-electron interactions, 
we assume a minimal model with a single interaction and 
1/Ti(r) cx l/r'^'^. 

The observed nuclear magnetization consists of nuclear 
spins with various configuration of randomly distributed 
impurities. The ensemble average is taken over all possi- 
ble impurity configuration. Using a probability distribu- 
tion, we obtain the nuclear magnetization 



{p{vi,t))Av = exp[-(^)host](ncxp[-; * 



exp[ 



'(rp )h0St] 
J-l 



]^{(l-c)+cexp[- 



/AV, 



exp[-( — )host 
J-i 



c 



d°r(l-exp[^])], (12) 



where c is the impurity concentration and the observed 
nuclear spins are assumed to lie on a space D dimension. 
The nuclei within a radius (wipeout radius) of r,, of the 
impurity are unobservable in conventional resonance fre- 
quency and time scale, so that they are excluded from 
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the integration. In Fig. (TTJ blue large circles indicate the 
wipeout regions of a radius r^. 

Taken into account a finite wipeout radius rc, the re- 
covery curve of eq. (jl2|) is expressed by 



pit) = p(0)exp[-( — )host 



where is wipeout number and <c is 7i(r = r^). After 
integrating eq. (|12p in the limit of — >0, the recovery 
curve is 

p[t) = p(0)exph(;^),ost -(-)"], 

= p{Q)eM-{^\ost ~ (14) 
Ji V n 

where an impurity-induced nuclear spin-lattice relaxation 
rate 

i_ oc ( ^ y/n helnhfS{S + 1) ^^^^ 

where n{—D/d) is 1/2 and r,„ is the decay rate of an 
impurity spin autocorrelation function. Thus, we obtain 
a theoretical recovery curve expressed by a product func- 
tion of exponential function times stretched exponential 
function, eg. p^ and eq. (fH|) . 

In ca. (fT4)) . p{0), (ri)host and ti are the fit parame- 
ters. The original recovery curve expressed by p{t) ~ 
p(0)exp[-V(Ti)host]n»[(l - c) + cexp(-t/ri(rO)] is a 
function of a lot of time constants of (Ti)host, Ti{ri), 
ri(r2), Ti(r3), • • ■ . For c <1, the product HJ- • ■] is ap- 
proximated by « exp[— c J dr{l — exp(— i/Ti(r)))], and 
then the spatial integral leads to the stretched exponen- 
tial function with the single time constant ti. Thus, 
only two time constants, (ri)host and ri, are obtained. 
(I'i)host is the Cu nuclear spin-lattice relaxation time due 
to the host Cu electron spin fluctuation via a hyperfine 
coupling. Tl is the impurity-induced nuclear spin-lattice 
relaxation time. The distribution of Ti is taken into con- 
sideration through Tl (r) but is not known a priori. Thus, 
we must assume the above-mentioned minimal form. 

These theoretical curv es are reproduced by more so- 
phisticated manner [l25j . In passing, one may find a 
mathem atica l ly sa me expression for diffusion-limitted re- 
laxation Il26l I127II . It should be noted that the diffusion- 



limitted relaxation is based on the nuclear spin-spin re- 
laxation in a long time t and the impurity-induced non- 
exponenital function in a short time t. Thus, the ap- 
plication of the form of eg. p^ itself is justified for any 
cases. 

Figure [12] shows the nonexponential ^^Cu(2) nuclear 
spin-lattice relaxation curves of Ni- and Zn-substituted 
Y1248 at T = 4.2 K and the least squares fits using 
eq. dni) (solid curves) [sij. In eq.([I31), p(0), (Ti)host, 



Ni vs Zn 
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FIG. 12: Right figures show nonexponential ®'^Cu(2) nu- 
clear spin-lattice relaxation curves of Ni (upper panel) and 
Zn (lower panel) substituted Y1248 at T = 4.2 K. The solid 
curves are least squares fits by eq. (|13|l based on impurity- 
induced NMR relaxation theory with wipeout effect. Left 
figures illustrate the Cu02 planes with the impurities and the 
wipeout regions. All the data are reproduced from [84[ |. 



tl, and iVc are the fit parameters. Both Ni and Zn im- 
purities induce nonexponential recovery curves. A differ- 
ence in Ni and Zn impurity effects is the wipeout number 
A'c and the size of the wipeout radius Tc- The wipeout 
radius around Zn is one order larger than that around 
Ni f84l.lT05l. 



B. SHRINK OF WIPEOUT REGION AROUN Zn 

Figure [13] shows the Zu substitution dependence of the 
wipeout number and the estim ated wipeout radius r^ 
for Y1237 and Y1248 at T 4.2 K [loj. The wipeout re- 
gion around Zn shrinks as Zn concentration is increased. 
This is associated with the Zn-induced change of host 
antiferromagnetic spin fluctuations. 

The size of the local magnetic state induced around Zn 
decreases as Zn concentration is increased, because the 
host antiferromagnetic correlation length ^af decreases 
with Zn (see below). On the other hand, the super- 
conducting virtual bound state (Andreev bound state) 
is induced around Zn and its size increases as Zn concen- 
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FIG. 14: Host <^3Cu(2) nuclear spin-lattice relaxation 
rate (l/rir)host and impurity-induced ^^Cu(2) nuclear spin- 
lattice relaxation rate 1/ri for optimally doped Y1237 with 
Ni i mpurity (a) and (b) [t^ and with Zn impurity (c) and 

(d) [ioi,[iQi • 



FIG. 13: (a)Wipeout number Nc and (b)wipeout radius Tc as 
functions of Zn concentration for Y1237 and Y1248 at T = 4.2 
K [l07l |. Wipeout region per Zn shrinks with Zn substitution 
at T = 4.2 K. 



tration is increase(i, because the upper critical field i?c2 
decreases with Zn and the supe rconducting coherence 
length ^sc increases with Zn 128| . A virtual bound stat e 
around Zn is actually observed by STS technique jl00l |. 
The core radius is estimated to be « ^sc- 

In the theory based on SO (5) symmetry, if ^sc is 
equal to ^af, the coexistence of superconductivity and 
antiferromagnetic long range ordering could be real- 
ized |l29l . Il30l |. In an underdoped regime, Zn recovers 
the coherency of magnetic correlation over the destruc- 
tive effect of magnetic dilution [I3l|] . Then, the effect of 
Zn on the size of the magnetic bound state may not be 
monotonic. Such recovery of the Nee l ord ering was also 
observed for La2-a;Sr2;Cu04 with Ni Il32l |. 



C. Zn VS Ni 

Although we believe that eq. ([13]) with a finite iVc repro- 
duces the actual recovery curve for impurity-substi tuted 
samples better than eg. p^ with a fixed iVc = jl05l |. 
there is a shortcoming of overestimation of above 
Tc [11]. Thus, we adopt the next best policy and ap- 
ply eq. pH) to obtain host and guest spin dynamics. 

Figure [TH shows (l/TiT)host and 1/ti using eq. (ITIl) 
for Ni- a nd Zn-substituted optimally doped Y1237 [Tg . 
[loi, [lo^. Figure [IS] shows (l/Tir)host and l/n us- 
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FIG. 15: 



Host 



^''Cu(2) nuclear spin-lattice relaxation 
rate (l/TiT)host and impurity-induced ®'^Cu(2) nuclear spin- 
lattice relaxation rate 1/ri for naturally underdoped Y1248 
with Ni imp urity (a) and (b) [s^ . Il04| and with Zn impurity 
(c) and (d) [lOallOg]. 



ing eg. (IT4l) f or Ni- and Zn-substituted underdoped 
Y1248 (84 ll0l[T06L[T09| . Both Ni and Zn induce non- 
exponential nuclear spin-lattice relaxation for the under- 
doped and the optimally doped samples. But, for low 
impurity concentration, the host spin fiuctuation is ro- 
bust for Ni substitution but fragile for Zn substitution. 
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The pseudo spin-gap behavior of 1 /Ti is robust both for 
Ni and Zn substitution. For high impurity concentra- 
tion, the strong wipeout effects on Cu NQR signals and 
the suppression of pseudo spin-gap are observed for both 
Ni- and Zn-substituted underdoped Y1248 [10^. This 
is consistent with the neutron scattering measurement 
for Zn- substituted oxygen deficient Y123 by Kakurai et 

ai [mt . 



D. Cu(2) NQR NEAR Zn 

Zn impurities induce low frequency satellite signals 
in the plane-site Cu(2 ) NQR spect ra for Zn -substituted 
Y1248 [lOi, [lOi [Hi and Y1237 [lOl, \lO^. From the 
relative intensity of Cu NQR spectra, the satellite sig- 
nals are assigned to Zn-neighbor Cu NQR lines. The Cu 
nuclear spin-lattice relaxation time of the satellite signal 
is shorter than that of the main signal in the supercon- 
ducting st ate o f Y1237 and Y1248 and even above Tc 
for Y1248 [loi|. Although Y1248 is a stoichiometric, ho- 
mogenous, underdoped electronic system, the Zn-induced 
inhomogeneous magnetic response in the Cu02 plane is 
more marked than that of the optimally doped Y1237. 
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FIG. 16: Zn-induce satellite signals and reduced main signals 
of ''^Cu NQR spectra for Y1237 and Y1248 reproduced from 
Itoh et al. [lOa] . 

Figure [TC] shows the Zn impurity substitution effects 
on 63Cu(2) NQR spectra at T = 100 K in the normal 
states of Y1237 (a) and Y1248 (b). Figure [T71 shows the 
63'65Cu(2) NQR spectra at T = 100 K for Y1237 (a) and 
for Y1248 (b) from the top to the bottom as functions 
of Zn concentration. The main signal intensity of the 
63Cu(2) NQR spectrum around 31.5 MHz for Y1237 and 
around 29.8 MHz for Y1248 rapidly decreases as Zn im- 
purity concentration is increased. The lower frequency 
tail and the broad satellite signal increases as Zn impu- 
rity concentration is increased. The solid and the dotted 
curves are numerical simulations using three Gaussian 
functions (Mi, M2, and S). By two Gaussian functions 
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FIG. 17; Zn-substitution effects on t he p lane-site Cu NQR 
spectra for Y1237 and Y1248 at 100 K [lo||. Individual NQR 
spectra are fit by three Gaussian functions (dashed curves) of 
Ml, M2, and S. Original integrated Cu NQR intensity of each 
pure sample is given by /(O) = Ml + M2 for Zn-free Y1237 
and Y1248. 



of Ml and M2, we can approximate the main spectra. 
The shaded Gaussian S is the satellite spectrum. 

In Fig. 1181 the relative intensity of the main spectrum 
Ima.iix{—Mi+M2) (upward solid triangles), the satel- 
lite spectrum /gate (=5) (open circles), and the total 
intensity(=Ml-|-M2-|-S') are plotted as functions of the 
plane-site Zn(2) concentrations for Y1237 (a) and Y1248 
(b). All the /main, /sato and the total intensity are nor- 
malized by /(O). In the upper panels, all Zn impurities 
are assumed to be substituted for the plane-site Cu(2); c 
= 3a;/2 for Y1237 and c = 2a; for Y1248. In the lower pan- 
els, the Zn impurities are assumed to be equally substi- 
tuted both for the chain-site Cu(l) and plane-site Cu(2); 
X for Y1237 and Y1248. 

The dashed curves are numerical calculations of Cu 
NQR intensity with various Zn configurations, using the 
binomial distribution function 



n(j)-k 



(16) 



where n{j) is the number of the jth nearest neighbor 
(nn) Cu sites. (1-c) is the probability that a plane site is 
occupied by a Cu ion but not a Zn ion, that is, the total 
Cu NQR intensity with Zn substitution. With respect 
to Zn configuration, the probability of finding Cu atoms 
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FIG. 19; Illustrations of (a) a '^^Cu(2) NQR spectrum and 
(b) a top view of the Cu02 plane with a Zn impurity for Zn- 
substituted Y1237. The shaded area in (a) is assigned to the 
Zn- neighbor Cu NQR signals. In (b), Zn locates at the center, 
open circles are the jth nn Cu sites near Zn, and solid circles 
are Cu sites away from Zn. 



FIG. 18: Zn-substitution dependence of integrated Cu NQR 
spectra using multiple Gaussian functions for Y1237 and 
Y1248 at 100 K. Red triangles are the main signal intensity. 
Blue circles are the satellite signal intensity. Open squares are 
the summation of the main and satellite signals. The dashed 
curves are numerical calculations using binominal distribution 
functions. The in-plane Zn concentration c is assumed to be 
c = 3x/2 for Y1237 and c = 2a; for Y1248 in upper panels, 
while in lower panels, c = a; is assumed both for Y1237 and 
Y1248. 



can be seen in the decomposition of 

1 = list ■ l2nd ■ Isid ' ' ' 

i=i 

j=l k=0 

The total number of probabiUty of finding Cu decreases 
as Zn concentration c is increased, that is (1-c). Keeping 
the lowest order terms of c(^ 1), we obtain 



(1-c) - (l-c)[]^„(,)Bfc(c) 

j = l k=0 

= (1 -c)[(l -c)^f=i"(^') 

N 



-f^n(i)c(l- 



0{c% (18) 



Figure [13 illustrates a ^^Cn{2) NQR spectrum (a) and 
the Cu02 plane with a Zn impurity (b) for Zn-substituted 



Y1237. The shaded area in (a) is the Zn-induced Cu 
NQR signal. In (b), Zn locates at the center, open circles 
are the jth nn Cu sites near Zn, and solid circles are 
Cu sites away from Zn. These are guides for the below 
calculations of NQR signal intensity. 

The calculated curves on the Cu NQR intensity are 
in what follows. In Figs. [THl the decreasing functions as 
increasing c are 



/(c)//(0) 



(1-C)4S0(C)-(1-C)5(1), 

(l-c)(4i?o(c))3-(l-c)i3(2), (19) 
(l-c)(4i?o(c))^ = (l-c)2i(3). 



These are the probability of finding the Cu atoms with 
all Cu atoms in their 1st nn positions (1), in their 12 
neighboring (up-to 3rd nn) positions (2), and in their 20 
neighboring (up-to 4th nn) positions (3). In other words, 
the diminished Cu NQR signals come from the Cu atoms 
with at least one Zn atom in their 1st nn positions (the 
1st nn wipeout effect) (1), in their up-to 3rd nn positions 
(wipeout effect up to the 3rd nn sites) (2), and in their 
up-to 4th nn positions (wipeout effect up to the 4th nn 
sites) (3). 

In Fig. 1181 the increasing functions as increasing c are 



/(c)//(0) 



(l-c)8Bi(c) 
(l-c)4Bi(c) 



8c(l 
4c(l 



H4), 



(20) 



These are the probability of finding the Cu atoms with 
one Zn atom in their 4th nn positions or in the 2nd and 
3rd nn positions (4), and in their j(=l, 2, 3)th nn posi- 
tions (5). 

In the upper panel of Fig. [THl the Zn-doping depen- 
dence of the main intensity for Y1237 (x <0.017) and 
Y1248 {x <0.022) is close to that of the case (3) within 
the experimental accuracy, where the Cu sites located up 
to the 3rd nn from Zn are unobservable (wipeout). The 
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Zn doping dependence of the total intensity for Y1237 
{x <0.017) and Y1248 {x <0.022) is close to that of the 
case (2), where the Cu sites located up to the 4th nn from 
Zn are unobservable (wipeout). For Y1237, the satellite 
signal is close to either of the Cu site at the 4th nn from 
Zn or the Cu sites at the 2nd and 3rd nn from Zn [the 
case (4)]. For x <0.017, the former assignment is consis- 
tent with those to the total and the main signals, whereas 
for a::=0.033, the latter is consistent with those to the to- 
tal and the main ones. For Y1248, the satellite signal is 
close to the Cu site at the 1st, 2nd, or 3rd nn from Zn [the 
case (5)]. However, this is inconsistent with the wipeout 
effect on the total intensity up to the 3rd nn. Thus, the 
satellite resonance of Y1248 may arise from the fraction 
of the Cu site at the 4th nn from Zn. 

The dif feren t assignment o f the satellite signal from 
the above [109| | is reported in [122| |. where the satellite is 
assigned to the 1st nn Cu site near Zn. The point that 
the satellite comes from Zn-neighbhor some Cu sites is 
consistent with each other. The wipeout effect on the 
Cu NQR main signals and the growth of the Zn-induced 
Cu satellite signals indicate that the satellite but not the 
main signal involves the Zn-neighbor Cu nuclei. However, 
another site assignment is proposed from theoretical cal- 
culation using cluster model, where the higher frequency 
side of the main Cu NQR signal involves the 1st nn Cu 
sites and the satellite signal is assigned to the 2nd nn 
Cu site to Zn [133]. No wipeout effect is supposed. The 
actual Cu NQR spectrum does not consist of sharp split 
lines but rather broad bands. That is a reason which 
makes it difficult to assign the signals. 



temperature is decreased. The difference in the tem- 
perature dependence of between the main and 
the satellite signals is marked even above Tc- The main 
^■^(l/Tir)host above about 100 K decreases with Zn sub- 
stitution. For low Zn concentration, the pseudo spin-gap 
temperature Ts is invariant, while for h igh Zn concen- 
tration, Ts rapidly decreases |l06l . Il09l |. The collapse 
of pseudo spin-gap suggests the occurrence of a quantum 
phase transition. In the underdoped regime, Zn impurity 
induces superconductor-insulator crossover or transition. 

The Cu(2) nuclear spin-lattice relaxation time probes 
the in-plane antiferromagnetic dynamica l spi n sus c epti- 
bility at an NMR/NQR frequency [HI, [Isi, [HE [l3^ 
and reflects homogeneity of the Cu02 pl ane. Rel ation 
between the Cu NQR relaxation results 'l06|, Il09j | and 
the ^^Y NMR results [uot is not straightforward. 

of the satellite signal is enhanced more than 
that of the main signal. This is an evidence for the Zn- 
induced virtual bound state JAj74_,94[via the locally en- 
hanced magnetic correlation [oj, [QoTIotI. [osj . The differ- 
ence in ^^(1/riT) of the satellite from the main indicates 
that Zn-induced "s taggere d moments" persist both above 
and below Tc [H, lOTl l98l. [9^ , where spatiall y extended 
potential scattering by Zn plays a key role |l37l . Il38| |. 
and that the pseudo spin-gap effect is suppressed near 
Zn [H^. 



V. NMR STUDIES OF VORTEX CORES 
A. REDFIELD PATTERN 



E. Cu(2) NUCLEAR SPIN-LATTICE 
RELAXATION TIME NEAR Zn 

The substitution of Zn impurity for Cu site broadens 
the Cu NQR spectrum and induces the low frequency 
satellite signals. Nuclear spin-lattice relaxation results 
from spin fluctuations through local nuclear sites in a 
real space [5l| . We are concerned about frequency distri- 
bution of Cu nuclear spin-lattice relaxation time and the 
temperature dependence. The detailed measurements of 
frequency dependence of Cu nuclear spin-lattice relax- 
ation curves revealed that the spin fluctuations depend 
on the distance from Zn impurity. 

Figure HOI shows and ^^(l/Tir)host of the 

satellite and main signals as functions of tempe rature 
for Zn-substituted Y1248, respectively [T06', ^109^. The 
relaxation rates were estimated by using ea. p^ . The 
main signals have both components of and 
^^(l/TiT)host, whereas the satellite signals have a single 
component of ^^(1/tiT ). Fo r Y1237, the satellite signals 
have both components [l09| . The difference in ^'^(1/tiT) 
between Y1237 and Y1248 is associated with the degree 
of underlying magnetic correlation and the size of the 
pseudo spin-gap. 

In the superconducting state, increases as 



Not only the impurity-induced magnetism but also 
magnetism inside vortex cores in the mixed states of high - 
Tc superc ondu ctors have attracted great interests Il29ll . 
Localized [l4Clj | and delocalized quasi-particle sates jl4l| 
are induced inside the vortex cores in conventional and 
unconventional superconductors, respectively. Vortex 
core magnetism and vortex core ch arge have been ex- 
plored from microscopic viewpoints jl42l . Il43l Il44| . By 
the STS techniques, scanning conductance spectra were 
measured around vortex cores and revealed the anoma- 
lous local density of states around the cores, indicat- 
ing the e xisten ce of internal electric and magnetic struc- 
tures [Tlfl45t. 

An NMR spectrum is known to show a characteristic 
pattern in the mixed sta tes o f type-II superconductors, 
that is Redfield pattern [14(| . This results from mag- 
netic field distribution of quantized magnetic flux pro- 
duced by local supercurrent against an external magnetic 
field. Quantized magnetic flux forms a triangular lattice 
or a square lattice. Each quantized magnetic flux has a 
vortex core. The magnetic flux spreads over a distance of 
a penetration depth Asci while the vortex core spreads 
over a distan ce o f the superconducting coherence length 

?sc(« Asc) [m. 

At the early stage, there had been a pessimistic view 
that NMR technique cannot detect the Redfield pattern 
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FIG. 20: Zn-substitution effects on ^'^Cu(2) nuclear spin-lattice relaxation rates at Cu NQR main and satellite signals ^"^(l/riT) 
and ®^(l/rir)host as functions of temperature for for Y1248. Open circles are ^^(l/TiT)host for pure Y1248. The dashed lines 
indicate the respective Tc's. Ts is the ps eudo spin-gap temperature defined by the maximum temperature of ''^(l/rir)host- AU 
the data are reproduced from [106| . Il09l ] 



for high-Tc superconductors, because of the long penetra- 
tion depth Asc- However, the problem is not whether Asc 
is long or short. In general, the width of Redfield pattern 
is at most the lower critic al field of Hd in a dilute vortex 
lattice (T3I Il48l Il49l . ll5Cll | . Sinceuon spin rotation mea- 
surements indicate Hci oc Tc [147j |. the Redfield pattern 
can be expected for higher Tc compounds. The problem 
is whether the linewidth i^nmr of NMR spectrum in a 
normal state is sharp eno ugh to detect additional broad- 
ening in the vortex state jl5l| . We have a criterion of 

Al/NMR < ^Hci, (21) 

ZTT 

where 7n is nuclear gyromagnetic ratio. 

Figure [21] illustrates internal magnetic field distri- 
bution h{r) due to quantized magnetic flux (a). Red- 
field pattern of NMR spectrum (b), and distribution of 
Knight shift Kspin{i^) due to d-wave local density of states 
(c) ^152]. 

A n external mag netic field is squeezed at the vortex 
core |l48l 1 1491 . Il50j | The squeezed magnetic flux density 
at the core produces a higher magnetic field than an ap- 
plied magnetic field, so that the shift at the nuclear site 
inside the core is larger than the shift in a normal state. 

Local density of states of electron s rec overs at the vor- 
tex cores of d-wave superconductors |l53l | . Then, not only 
the distribution of an applied magnetic field but also the 
distribution of Knight shift due to anisotro pic lo cal den- 
sity of states contributes the NNR pattern [l52j | . 

Using a vortex core shift '"^ i^c, field (V) >0 due to the 
enhanced field and a negative shift ^''-ftrc,dia(C)<0 due to 
the diamagnetic field, we obtain the resonance frequency 
j'(V) inside and i'{G) outside the vortex cores, 



i/(V) = j.o[l ifc,spi„(V)][l ifc,flcid(V)], (22) 

and 

KC) = t.0[l+^^^c.dia(C)]. (23) 

One should take into consideration the combined effect 
of magnetic flux squeezi ng a nd of the finite positive spin 
shift at the vortex core [152!|. 

The Redfield pattern was actual ly ob served in type - 
II supercon ductors, e. g. vanadium 'l46'|, LiTi204 [l54| . 
and MgB2 [l55l . Il56l |. In the high-Tc cuprate super- 
conductors, the Redfield patterns were also observed, 
e.g. by I'^O NMR of optimally doped Y1237 and un- 
derdoped Y1248 [US, llH, EHI, HB^, EE^ and ^ostI 
NMR of overdoped TlsBazCuOg+a (T12201) (T^ = 85 
K) [m [m, [l6|. The frequency distribution of ^^Ti 
in the Redfield pattern was o bserved for the first time 
by Curro et al for Y1237 [l57j |. This is the beginning of 
site-selective NMR studies of vortex core magnetism. 

In passing, the evidence of static antiferromagnetic 
vortex cores at Tn = 20 K and at = 2 T was re- 
ported for T12201 |l63l |. But, the irreversible magneti- 
zation line was identified at Ti„ = 20 K and at H = 
2 T, and the slowing down effect of vorte x me lting or 
solidification was observed in 209^1 NMR |l52l |. From 
magnetic field dependence of ^°^Ti, a direct process of 
overdamped mo tion of pancake vortices was reported to 
play a key role [l64l |. Whether the static antiferromag- 
netic ordering takes place inside vortex cores, one needs 
further studies. 
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FIG. 22: From top to bottom: zero field '''^Cu (nuclear spin 
I = 3/2) NQR spectrum at T = 100 K, Zeeman-split "^^Cu 
NQR spectra at T = 100 K and then 10 K at an e xternal 
magnetic field of -ff = 96 mT along the c-axis [l5l| . The 
Zeeman-split directions of '''^Cu NQR lines (solid lines with 
UQi and i^Qi) are denoted by arrows. The diamagnetic shift 
is superimposed on the original broad NQR lines. 



FIG. 21: Illustrations of magnetic field distribution (a), Red- 
field pattern of NMR spectrum (b) , and dist ribu tion of Knight 
shift due to d^2_y2- wave vortex cores (c) [152( |. Red shaded 
area indicates "Redfield pattern" of distributed Knight shift. 
One should note that magnetic field distribution due to su- 
percurrent is at most Hd . 

As to the plane-site Cu(2) NMR , the ideal condition of 
eq.dH]) is not satisfied for Y1237 [165| |. However, in the 
attempts to detect the spatial distribution of antiferro- 
magnetic spin correlation, we obtain a few results of Cu 
NMR studies of the vortex cores. 



100 K and 10 K (< Tc)- The magnetic field of = 96 mT 
was applied along the c axis, the maximum principle axis 
of electric field gradients. An intervortex distance at 96 
mT is about 1800 A, longer t han t he penetration depth 
\ab{T -^0 K) = 1100-1300 1 [mI. Hal is estimated to 
be about 110 mT at T ^0 K [iTlj . 

Magnetic field splits the ^'^Cu NQR spectrum with a 
peak frequency i^nqr into two li nes with lo wer and higher 
frequency peaks {i^qi and i'q2) |l72l Il7d |. By definition 
of 63-y* = (i-^ii:^j63^^ (63^^ ^ 11 285 MHz/T and the c- 

axis Knight shift Kcc), the transition frequencies of z/qi 
(/z=3/2^1/2) and iyQ2 (/z=— 3/2^— 1/2) are given by 



B. ZEEMAN-PERTURBED Cu NQR 

Let us present Zeeman-perturbed Cu NQR studies of 
vortex cores for optimally doped Y1237 at a low magnetic 
field of = 96 mT ^151). The plane-site Cu(2) NQR 
spectrum is known to be one of the sharpest among the 
reported NQR spectra of other cuprate superconductors. 
The Cu(2) nuc lei are cou pled with nearly uniaxial electric 
field gradients jl66l . At high magnetic fields, various 
vortex pinning effects are observed from measurements of 
hysteresis curves of magnetization [168]. At low magnetic 
fields near Hd, however, a large diamagnetic response is 
observed in any case. The vortex lattice at a low mag- 
netic field is in a quasi-long range ordering state [l69|. 
In the Cu NMR experiment for Y1237, the low magnetic 
field yields only a negligibly small Cu Knight shift [l70j | 
but is expected to cause a large diamagnetic shift. 

Figurel^shows a zero- field ^'^Cu NQR power spectrum 
at 100 K (> Tc) and the Zeeman-split NQR spectra at 



VC12 — J^NQR 



63 



(24) 



One should note that Kcc involves both the Knight shift 
and the superconducting diamagnetic shift. From vqi 
and vq2i one can estimate the quadrupole frequency 
v^QKiH) at a magnetic field H and the local magnetic 
field Sh[={l+Kcc)H] by 



'^NQR(-ff) - 

6h ^ {vQi 



{VQI + ;^Q2)/2, 

-J^Q2)/2'^3^„. 



(25) 



The Cu NQR spectrum at 100 K is broadened by the 
distribution of electric field gradients. The split ^^C\x 
NQR lines is broadened at 10 K larger than 100 K. The 
local field due to the Knight shift jl70l | is negligible at 
the magnetic field of iJ = 96 mT. Thus, the additional 
broadening at 10 K must be due to the superconduct- 
ing diamagnetic shift ^'^Kc^dia. and the distribution of the 
electric field gradient vnqk{H) in a mixed state. 
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FIG. 23; (a) Frequency distribution of Cu nuclear spin- 
lattice rate divided by temperature l/^^TiT at T = 10 and 5 
K (left axis) in the Z eema n-perturbed ®^Cu(2) NQR spectrum 
of VQ2 (right axis) [l5l| . The arrow indicates a minimum 
of Xj^^TxT for nuclei close to the saddle points in a vortex 
lattice, (b) Temperature dependence of the ®^Cu nuclear spin- 
lattice relaxation rate divided by temperature l/^'^TiT across 
the split NQR spectrum 151]. The arrow indicates a shallow 
minimum of l/^'^TiT at 10 K for the 29.5 MHz nuclei close 
to the vortex cores. 



From the estimated 31.54 < i/NQR(96mT) <31.58 MHz 
and ;/nqr(0) = 31.54 MHz at 10 K, we obtain the upper 
limit of the field-induced quadrupole frequency 

< KQR(96mT) - i/nqr(0)] < 40kHz. (26) 



Cu directly probes the in-plane antiferromagnetic corre- 
lation. The frequency dependence of 1/^^Ti is primar- 
ily explained by Doppler shift in q uasipartic le excitation 
spectrum around the vortex cores [l53l . Il75j . The super- 
current around a vortex core induces Dop pler shift in the 
local quasiparticle energy spectrum [l7g. Thus, the fre- 
quency dependence of 1 j^'^Ti , similar to that of the high 
field 1/^^Ti, indicates the importance of the Doppler shift 
in quasiparticle energy spectrum and the absence of the 
static antiferromagnetic vortex cores at 96 mT. 

Figure I^HTb) shows l/^^TiT as a function of temper- 
ature across the split NQR spectrum of h'Q2 and at zero 
field NQR. The signals at 30.4-30.6 MHz and at 29.5 
MHz in the split NQR spectrum come from the nuclei 
close to the saddle points and close to the vortex cores, 
respectively. The dashed line indicates l/TiT oc of d- 
wave pairing. The temperature dependence of 
at each site can be understood primarily by the spatial 
dependence of the local density of states with the Dopple 
shift ^EM- The minimum of l/^^^iT at 10 K for nuclei at 
29.5 MHz is not understood simply by the local density 
of states. 

Two theoretical explanations have been proposed. One 
is the competition mechanism in the scattering of quasi- 
particle and antiferromagnetic spin fluctuations [l75|. 
The other is the mismatch effect of angular dependent 
local density of states around a vor tex core and the ge- 
ometry of magnetic flux lattice 

In the presence of a supercurrent, the local antiferro- 
magnetic spin fluctuations have two fold contributions to 
l/^^TiT: one is a scattering of a Bogoliubov particle at 
a spin-fermion vertex, and the other is simultaneous cre- 
ation or annihilation of two Bogoliubov particles 17{|. 
Competition of two processes yields minimum behavior 
as functions of temperature and frequency. 

The periodicity of vortex lattice, i.e. triangular or 
square lattice, depends on the magnitude of an external 
magnetic field H. Mismatching between the magnetic 
flux shape and the local density of states of the quasi- 
particles yields the c hara cteristic dependence of temper- 
ature and frequency [l77j . 



The low magnetic field of 96 mT does not enhance the 
vortex charge, compared with that of 9.4 T 174| . There 
is no robust evidence of vortex charging effect on NQR 
frequency beyond the experimental accuracy. 

Figure l^a) shows the frequency dependences of the 
^^Cu nuclear spin-lattice relaxation rate divided by tem- 
perature if^^TiT at T = 10 and 5 K in the Zeeman- 
perturbed ^^Cu(2) NQR spectrum of vq2 (right axis). 
Above Tc, shows no appreciable frequency depen- 

dence, whereas below Tc, it shows a strong dependence. 
This is an evidence for the site-selective measurement of 
^^Ti in the Zeeman-split NQR spectrum and for the less 
effect of nuclear spin diffuison. 

The frequency dependence of 1 /^^Ti is similar to those 
of the plane-site ^^0(2, 3) and the apica l "'^^0(4) at high 
magnetic fields of = 9-37 T [lHlllOl. The plane-site 



VI. CHEMICAL PRESSURE, PHYSICAL 
PRESSURE, AND SITE DISORDER 

It has been believed that the site disorder out of Cu02 
planes does not cause serious damage to the electronic 
states of the Cu02 planes. However, the reason why 
chemical pressure effect on Tc and the electronic states 
by element substitution is not the same as the physical 
pressure one has been often ascribed to the disorder ef- 
fect. The relation among site disorder, chemical pressure 
and physical pressure effects has been poorly underdtood. 

Y( Bai .,:S.7,)2Cu408 with x < 0.40 can be synthe- 
sized (l79l |. For Y1248, one can expect two effects of Sr 
substitution for Ba sites without change of oxygen con- 
tent; one is chemical pressure and the other is crystalline 
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potential disorder. The size of Sr^+ ion is smaller than 
that of Ba^+ by about 10 %, so that the substituted Sr 
ions make local strains and decrease the volume. Th e lat - 
tice constants actually shrink with Sr substitution jl79l |. 
similarly to physical pressure effect [180]. Thus, the Sr 
substitution introduces crystalline potential disorder and 
chemical pressure to the lattice. 

The effect of hydrostatic physical pressure o n the spin 
dyn amics is similar to the carrier doping effect |l8ll . Il82l 
Il83j |: The physical pressure effects on and on the 
pseudo spin-gap temperature Tg agree with the carrier 
doping effects for Y1248, Tc increases fisi. [l85l bu t 
decreases |l83j |. as well as with doping Ca |18lTl86l |. 

For Y12 48, th e Sr substitution for Ba site scarcely in- 
creases Tc [iZSl, although the unit cell volume shrinks. 
Thus, one can suspect randomness effect on the actual Tc 
as well as nonmagnetic impurity Zn substitution effect. 
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FIG. 24: Zero field '^^•'^''Cu NQR frequency spectra of the 
chain-site Cu(l) (left side) and the plane-site Cu(2) (right 
side) for Y(Bai_^Sr:r)2Cu408 of cc = 0.00 and 0.40 at an 
ambient pressure of P = 0.1 MPa [l78j . The red shaded Cu(l, 
2) NQR spectra were measur ed u nder hydrostatic physical 
pressure of P = 2 and 3 GPa [iS^ . 

Figure [Ml shows ^'^'^^Cu NQR spectra of the plane-site 
Cu(2) (a) and the chain-site Cu(l) (b) for x = 0.00 and 
0.40 at T = 4.2 K 178]. For comparison, the red shaded 
63,65(^^ NQR spectra for a pure Y1248 under hydrostatic 
physica l pre ssure of P — 2.0 and 3.0 GPa are also pre- 
sented |l82| . By Sr substitution for Ba site, the Cu(2) 



NQR spectra are shifted to higher frequencies, whereas 
the Cu(l) NQR spectra are shifted to lower frequencies. 
Since the directions of these s hifts a re in parallel to those 
under the physical pressure jl82l . Il87j ]. However, one 
should note that the degree of shifts of the Cu(2) NQR 
spectra is different from that of Cu(l) betwe en the Sr - 
substitution and the physical pressure effects |l82l . Il87l ] . 
This difference indicates that local compression due to 
the internal pressure of Sr is different from that due to 
the physical pressure. 

Both the line widths of Cu(l) and Cu(2) spectra are 
broadened by Sr substitution and then the disorder of 
the crystalline potential due to Sr in BaO layer. For such 
a broadened spectrum, one would expect an inhomoge- 
neous local density of electron states and an inhomoge- 
neous electron spin dynamics. However, no impuri ty ef - 
fcct on nuclear spin-lattice relaxation was observed [l78j ] . 
Both Cu(l) and Cu(2) nuclear spin-lattice relaxation 
curves were single exponential functions. Also, no Curi e 
term in uniform spin susceptibility was observed [178| ]. 
Theses results are sharply in contrast to the effects of 
in-plane impurity substitution. It remains to be a mys- 
tery why the chemical pressure of Sr substitution does 
not increase Tc. 



VII. CONCLUSION 

Local electronic states near the in-plane impurities and 
inside the vortex cores have been studied by Cu NQR 
and NMR techniques for high-Tc superconductors Y123y 
and Y1248. The characteristic NQR and NMR spectra 
and the frequency distribution of Cu nuclear spin-lattice 
relaxation time provide us with rich information of mag- 
netic correlation near Zn impurity and inside the vortex 
cores. 
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